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ABSTRACT 

Context. Energetic pulsars power winds of relativistic leptons which produce photon nebulae (so-called pulsar wind 
nebulae, PWNe) detectable across the electromagnetic spectrum up to energies of several TeV. The spectral energy 
distribution has a double-humped structure: the first hump lies in the X-ray regime, the second in the 7-ray range. 
The X-ray emission is generally understood as synchrotron radiation by highly energetic electrons, the 7-ray emission 
as Inverse Compton scattering of energetic electrons with ambient photon fields. The evolution of the spectral energy 
distribution is influenced by the time-dependent spin-down of the pulsar and the decrease of the magnetic field strength 
with time. Thus, the present spectral appearance of a PWN depends on the age of the pulsar: while young PWNe are 
bright in X-rays and 7-rays, the X-ray emission of evolved PWNe is suppressed. Hence, evolved pulsar wind nebulae may 
offer an explanation of the nature of some of the unidentified VHE 7-ray sources not yet associated with a counterpart 
in other high-energy ranges. 

Aims. The purpose of this work is to develop a model which allows to calculate the expected X-ray fluxes of unidentified 
VHE 7-ray sources considered to be PWN candidates. Such an estimate helps to evaluate the prospects of detecting 
the X-ray signal in deep observations with current X-ray observatories in future studies. 

Methods. We present a time-dependent leptonic model which predicts the broad-band emission of a PWN according to 
the characteristics of its pulsar. The values of the free parameters of the model are determined by a fit to observational 
VHE 7-ray data. For a sample of representative PWNe, the resulting model predictions in the X-ray and 7-ray range 
are compared to observations. 

Results. The comparison shows that the high-energy emission of identified PWNe from different states of evolution is 
predicted correctly by the model. 

Key words. Pulsars: general; Gamma rays: general; ISM: supernova remnants; ISM: individual objects: (MSH 15—5,2; 
HESS J1825-137; HESS J1837-069); Radiation mechanisms: non-thermal; 



1. Introduction 

During the last decade, a new generation of Imaging 
Atmospheric Chercnkov Telescopes (IACTs) has discovered 
an increasing number of Galactic sources emitting VHE 
(very high energy, E > 100 GeV) 7-rays. Among these, pul- 
sar wind nebulae (PWNe) form the most abundant class. 
Such nebulae arc usually associated with the non-thermal 
emission from a magnetized plasma of relativistic particles 
fed by an energetic pulsar. In current models, the plasma 
is thought to consist m ainly of energetic leptons (see e.g. 
iGaensler fc Sl anc. 200(|) which emit non-thermal radiation 
over a wide energy range. Interacting with magnetic fields, 
the leptons produce synchrotron radiation up to several 
kcV. In addition, low-energy photons, e.g. from the cos- 
mic microwave background (CMB), can be up-scattered by 
the energetic leptons to very high energies via the Inverse 
Compton effect. Therefore, the emission in X-rays and 
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VHE 7-rays is tightly li nked, emerging from the same lep- 
ton population (see, e.g. lGelfand et all I2009T ). 
The second largest population of Galactic VHE 7-ray 
sources consists of unidentified sources without an un- 
ambigious counterpart at other wavelengths (see, e.g., 
lAharonian et all I2008D . However, in many cases an ener- 
getic pulsar can be found in the vicinity, suggesting a possi- 
ble connection between these unidentified objects and pul- 
sar wind nebulae. Provided that the rotational period and 
its first time derivative can be measured, e.g. by radio ob- 
servations, the spin-down energy loss of the pulsar can be 
estimated and hence the viability of the pulsar as an en- 
ergy source of the nebula can be investigated. In addition, 
in the PWN scenario of broad-band emission by energetic 
leptons a spatial association of the VHE 7-ray source with 
an X-ray nebula counterpart is expected. There are mainly 
two issues complicating this identification scheme: In some 
cases, PWNe are slightly displaced from the pulsar, which 
may result from an inter action with the supern ova remnant 
reverse shock (see, e.g., iBlondin et all |2001[ ) and from a 
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prop er motion of the pulsa r , gain ed from a kick at its birth 
( e.g. Ivan der Swaluw et aD . l2004[ ). Furthermore, in particu- 
lar for older systems the X-ray emission becomes harder to 
detect since the energetic leptons injected in earlier epochs 
have been cooled and, at the same time, the supply of fresh 
leptons is reduced. Moreover, the synchrotron emission by 
the freshly injected leptons is suppressed because the mag- 
netic field strength decreases with time. Since the accumu- 
lated less energetic leptons can still produce VHE 7-rays 
via Inverse Compton (IC ) scattering, such evolved PWNe 
have been proposed (e.g. Ide Jager fc Diannati-Ata'i I2009D 
as an explanation of some of the as yet unidentified VHE 
7-ray sources. 

In this work we introduce a time-dependent leptonic model 
of the non-thermal emission of PWNe (Section^ and ap- 
ply the model to PWNe of different evolutionary states 
(Section[3]). For each individual source, the free parame- 
ters are fixed by fitting the model to the VHE 7-ray data. 
Subsequently, we show that the X-ray emission of these 
objects is predicted correctly by the fitted model. Hence, 
in future studies it may serve as a means to estimate the 
X-ray flux of unidentified VHE 7-ray sources in a PWN 
scenario, allowing to evaluate the prospects of detection in 
deep observations with current X-ray observatories. 

2. The Model 

In this Section we introduce a leptonic model describing the 
time evolution of the non-thermal radiation from PWNe. 
The time dependence of the energy output E of the pulsar, 
which derives from the slow-down of the rotation, has to be 
taken into account: 



E 



dE lx 



dt 



(1) 



Following iPacini fc Salvatil (119731 ). the energy output 
evolves with time as 



(2) 



E(t) =£ li- 
re 



where Eq = E(t = 0), tq denotes the spin-down timescalc 
of the pulsar and n the braking index. The latter has been 
measured only for a few young pulsars (|Camilo et all . 120001 
and references therein). Such a measurement exists, for in- 
stance for PSR B1509— 58, which is one of the sample pul- 
sars discussed in Section [3J Thus, for the modeling of the 
PWN associated with P SR B1509-58 we used the mea- 
sured value of n = 2.839 dLivingstone et aL . 2005tl. For th e 
other cases we adopted n = 3 ((Manchester &: Tavlorl . |1977|) , 
corresponding to spin-down via magnetic dipole radiation. 
The spin-down timescale tq is defined by 



To 



2r c 
n- 1 



with Pq and P being the initial and the current period, 
respectively, and r c = Pj (2P) the characteristic age of the 
pulsar (P denoting the time derivative of the rotational 
period). For n = 3 and Pq <C P the present true age T of 
the pulsar corresponds to the characteristic age r c , whereas 
for other cases it can be calculated as 



T = 



P 



(n- 1)P 



1 - 



(4) 



P and P can usually be derived from radio observations, 
while the initial period will be treated as a free parameter 
of our model. 

In the following, the evolution of the non-thermal emis- 
sion of a PWN is calculated in discrete time steps with an 
adaptive step size of St. In each time step only a fractional 
amount AE p of the energy output of the pulsar is con- 
verted into relativistic leptons, i.e. electrons and positrons. 
Assuming the corresponding conversion efficiency 77 to be 
constant over time, AE p (t) is determined as 



AE p (t) = T) 



t+st 



E(t')dt' 



(5) 



for a time interval [t, t + St] and with 77 G [0; 1]. The conver- 
sion efficiency is strongly correlated with Pq , such that the 
quality of the fit does not benefit from an additional free 
parameter. Therefore 77 is fixed to the value of 0.3, which is 
in agreement with e.g. the modelin g re sults for MSH 15— 5 2 
carried out bv lSchock et al.l (j2010ft and lZhang et all ([20081 ) . 
We assume that the differential energy spectrum of the in- 
jected leptons can be described by a simple power law 
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Assuming that the spectral shape does not change within a 
time bin, $o(t)-, denoting the normalization of the distribu- 
tion at 1 TcV, can be calculated by integrating the injection 
spectrum over energy: 



AE p (t) 



dJVi, 



dE 



-(E,t) dE. 



(7) 



We only consider leptons injected into the pulsar wind in 
the range between E min = 0.1 TeV and E max = 1000 TeV, 
well suited to accomodate the VHE 7-rays as well as the 
X-ray emission from PWNe. Leptons with energies outside 
this range do not significantly contribute to the emission 
in the considered photon wavebands. Given a differential 
yield of leptons dN(E, t — St)/dE with energy E at a time 
t — St the number of leptons remai ning after cooling a t the 
time t can be calculated. Following I Zhang et al.l (120081 ). the 



cooling of the lepton population during a time step St is 
implemented in the model by an exponential function: 



dN r , 



Died 



dE 



(E,t) 



dN 
dE 



(E, t-St)- exp 



St 



T eS (E,t) 



(8) 



This approach uses an effective cooling time scale 
T esc taking into account both synchrotron and 



'off 



'syn 



escape losses. The res pective time scales T syn and r GSC are 
likewise adopted from lZhang et all (|2008h : 



T syn {E,t) = 12.5- 



B(t) 



10 fiG 



E 



r esc (E,t) =34- 



B(t) 
10 fxG 



E 
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kyr, (10) 



where R(t) and B(t) describe the time evolution of the 
PWN radius and the magnetic field strength inside the 
PW N, respectively. For evo lved PWNe, R(t) is given by 
(see iGaensler fc Sland . [20061 and references therein): 



R(t) 



b ■ t 3 / 10 
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for 
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where the coefficients a and b can be calculated using the 
present-day size of the PWN. Likewise, we assume that 
the radial extent of each freshly injected lepton population 
can be described by Eq. [TTJ where the age of the PWN is 
replaced by the corresponding age of the leptons. This ap- 
proximation leads to an onion-like structure of the PWN 
model. 

The magne t ic fie ld strength B(t) is adapted from 
IZhang et all (|2008ft : 



B(t) 



B 



ISM 



(12) 



-Bism represents a time-independent component of 3 fiG 
to account for the magnetic field strength of the ambient 
medium. The parameter a is fixe d to a value o f 0.5, w hich 
was also found in the studies of I Zhang et all (120081 ) and 



iQiao et all (pOOl for five PWNe (among them MSH 15-52 
and HESS J1825— 137 which are investigated in this work). 
Finally, Bq, the initial magnetic field strength inside the 
PWN, is a free parameter. All in all, the model has two 
free parameters (Pq and Bq) defining the starting condi- 
tions of the PWN evolution. 

Having established the framework for cooling and injection 
processes, we can calculate the number of leptons with en- 
ergy E present in the nebula at a time t + St. This number 
comprises leptons injected and cooled until time t as well 
as freshly injected leptons between t and t + St: 



dE y ' 



St) 



dN, 



cooled 



dE 



(E,t) 



dW 



dE 



(E,t + St). (13) 



By iteratively evaluating Eq. Q21 it is possible to de- 
termine the energy distribution of the leptons inside the 
PWN at an arbitrary time. Based on this distribution, the 
corresponding photon population can be calculated, with 
synchrotron radiation and IC scattering as the most rele- 
vant emission processes in the considered energy rangcQ- 
A detailed account of t hese mechanisms can be found 
in IBlumenthal fe Gouldl (| 19701) . The target photon fields 
considered for IC scattering - CMB, starlight and in- 
frared photons - are adopted from the GALPROP code 
([Porter fc Strong! . I2005D . 

As a first step, we can use the model to study the develop- 
ment of the spectral energy distribution (SED) with pro- 
gressing age for a generic PWN system. The SED shown 
in Fig. [T] is based on the characteristics of the pulsar 
PSRJ1826-1334 and its nebula (see Table [J), represent- 
ing an example of an evolved PWN. The free parameters 
are exemplarily set to P = 30 ms and B = 50 ^G. 
Since the magnetic field strength decays strongly with time, 
the X-ray emission is suppressed for high PWN ages. At the 
same time energy-dependent cooling effects become visible 
in the 7-ray band, reducing the emissions in the VHE 7-ray 
range and shifting the peak to lower energies. 
As a second application of the model, we investigate the 
contribution of leptons from different injection epochs to 
the current photon SED, as shown in Fig. [5]for our generic 
PWN. Summing up the SEDs from leptons of all epochs 
results in the emission visible today. The modeling shows 



1 We neglect a synchrotron self-Compton (SSC) scattering 
component in the VHE 7-ray spectrum, since this work is fo- 
cused on evolved pulsars, whereas SSC is relevant mostly for 
the highly magnetized PW Ne of very young a nd energetic pul- 
sars, e.g. the Crab Nebula ijMever et al.l . l2010f ). 




1 (T 1 10 1 
Energy [TeV] 

Fig. 1. Typical evolution of the modeled spectral energy 
distribution of a PWN with time. The color scale repre- 
sents the age of the PWN, starting with a young system 
(500 years, yellow) and proceeding in equidistant steps on 
a logarithmic time scale to an old system (200 kyr, dark 
red). 
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Fig. 2. Photon SED (black broken line) of a generic middle- 
aged (approximately 20 kyr) PWN decomposed into con- 
tributions by leptons from different injection epochs (solid 
colored lines) . The same parameters as in Fig. [T] were used 
for this example. The grey areas represent the energy ranges 
covered by current X-ray and VHE 7-ray observatories. 



that mainly the youngest leptons account for the X-ray 
emission while we observe the injection history of the pul- 
sar in VHE 7-rays. This can be understood as the en- 
ergy of the leptons causing the synchrotron radiation in 
the considered X-ray regime is higher than the energy of 
the lep tons causing the VHE 7-ray em ission via IC scat- 
tering (|de Jager fc Diannati-Atal . I2009D . Numerous highly 
energetic leptons are present in young populations. For 
older populations the number of highly energetic leptons 
has been reduced significantly by cooling processes (mainly 
synchrotron radi ation), which are especially e fficient at high 
energies (see e.g. IBlumenthal fc Gouldl . 1 1970f) . Thus, for the 
older populations the synchrotron peak is shifted to lower 
energies, out of the observational X-ray range. 
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Table 1. Overview of the selected PWNe and their associated pulsars. The list is sorted by increasing characteristic age, 
representing different evolutionary states. The properties of the pulsars (characteristic age r c , current period P, current 
spin-down luminosity E and distance d) are taken from the ATNF pulsar database* (Manches ter et al.l . l2005| ). References 
for the H.E.S.S. sources: ^Aharonian et all (|2005h . [2 lAharonian et all (|2006af ). t 3 lAharonian et all (|2006bh . 



VHE Source 


equiv. VHE source radius 
[arc min] 


Pulsar 


[kyr] 


P 

[ids] 


[erj 


E 


d 

[kpc] 


MSH 15-52 


i] 


11.5 


PSR B1509 


-58 


1.55 


151 


1.8 


■10 37 


5.81 


HESS J1825- 


-137 [2] 


44.0 


PSRJ1826- 


1334 


21.4 


101 


2.8 


■ 10 36 


4.12 


HESS J1837- 


-069 131 


13.7 


PSRJ1837- 


0604 


33.8 


96 


2.0 


■ 10 36 


6.19 



*URL: http : //www . atnf . csiro . au/research/pulsar/psrcat/ 
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Energy [TeV] Energy [TeV] 



(a) MSH 15-5,8 (b) HESS J1825-137 




Energy [TeV] 

(c) HESSJ1837-069 

Fig. 3. Spectral energy distributions for the three sources listed in Table [TJ The black lines show the modeled SEDs 
resulting from a fit to the VHE 7-ray data. The violet bands depict the uncertainty of the model based on errors on the 
fit parameters and their correlations. The blue dotted bands denote the model prediction of the X-ray emission calculated 
for the published analysis regions (compare Table 2]). H.E.S.S. data (red filled circles) are presented along with their 1 a 
statistical errors. Red filled bands show the X-ray data. Where available we included 7-ray data (orange) from the Fermi 
Large Area Telescope (Fermi- L AT), as well. For better visualization we show higher- level error bands for the X-ray data 
and for the modeled SEDs. References for the VHE 7-ray and X-ray dat a can be found in Tables Q] and [ H respectively . 
Fermi data for MSH 15-52 and for HESS J1825— 137 are adopted from lAbdo et~aTl (f2010h and l&ondinet~al] (|20ll . 
respectively. 



3. Applications 

Based on the considerations presented above, evolved 
PWNe offer a potential explanation for a significant frac- 
tion of TeV sources which have remained unidentified up to 
now. In this scenario, an old pulsar may be surrounded by 



a relic TeV PWN detectable with current IACTs. However, 
due to the low present energy output of the pulsar the neb- 
ula does not contain enough high-energy leptons to pro- 
duce a strong X-ray counterpart to the TeV PWN (com- 
pare Figs. [1] and [2]). This view is in accordance with de- 
tailed studies presented, e.g., bv lde Jager fc Diannati-Ataj 
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Table 2. Application of the model to the selected PWNe. The model was fit to the VHE 7-ray data and the fit results 
of the free parameters Pq and Bq are listed. The interval in parentheses denotes the parameter boundaries during the 
optimization procedure. 



Source 


Po 
(5- 


[msj 
-P) 


Bo 
(5- 


200) 


X 2 /n.d.f 


MSH 15-52 




43.5 


± 2.5 


85.0 


±8.3 


15.1/12 


HESS J1825- 


-137 


28.4 


±0.8 


26.9 


± 1.5 


30.5/9 


HESS J1837- 


-069 


39.6 


±3.7 


14.4 


±3.0 


10.0/10 



Table 3. Comparison of predicted and measured VHE 7-ray flux for the modeled PWNe. References can be found in 
the caption of Table [T] 



Source 


Energy threshold 
[TeV] 


Flux above threshold [10~ 12 cm -1 s _i ] 
Model Prediction Measured 


MSH 15-52 
HESS J1825- 
HESS J1837- 


137 
069 


> 0.28 

> 0.27 

> 0.2 


19.8 ±1.1 22. 5 a 
73.8 ± 1.8 82.8 ±2.2 
25.2 ± 2.6 30.4 ±1.6 



a No statistical errors provided. 



Table 4. Comparison of predicted and measured X-ray emission for the modeled PWNe. For each X-ray analysis the 
published energy range and the radius of the used analysis region are listed. In the case of MSH 15— 5,2, where the authors 
provide several ring-shaped analysis regions, we co mbined the innermos t t hree rings in order to use c omparable extraction 
areas for all sources. References to the X-ray data: ^ ISchock et al.l (|2010t ). [ 2 jUchivama et all (|2009l ). t 3 lGotthelf Sz Halpernl 

posh . 



Source 


Analysis region 
Radius [arcmin] 


Energy range 
[keV] 


Px 

Model Prediction 


Measured 


Index" 
Model Prediction 


Measured 


MSH 15-52 
HESS J1825- 
HESS J1837- 


i] 

-I37P] 
-069 131 


0.5-2.3 
1.5 
1.0 


0.5 - 9 
0.8 - 10 
2-10 


22.1 ± 2.0 
1.82 ±0.07 
2.5 ±0.4 


24.9 ±0.7 
1.71 ±0.47 
1.0 d 


1.54 ±0.04 
1.53 ±0.01 
1.57 ±0.13 


1.66 ± 0.02 c 
1.69 ±0.09 
1.6 ±0.4 



a Energy fluxes are given in units of 10 -12 ergcm -2 s _1 . 
b Photon index of a power-law model. 

c Photon index from the innermost ring of the X-ray analysis which exhibits the highest energy flux. 
d No statistical errors provided. 



(|2009t) and lMattana et all (|2009t ). For some of the uniden- 
tified VHE 7-ray sources where an evolved PWN scenario 
appears likely, deep observations performed with current 
satellite observatories may yet reveal X-ray counterparts 
despite the relative weakness of the expected X-ray emis- 
sion. The model presented in this work allows to select 
suitable candidates based on an estimate of the required 
exposure for a detection in the X-ray regime. In order to 
investigate the reliability of the model we applied it to three 
selected PWNe, which are listed in Table [T] The motivation 
for this selection was to sample PWNe from different states 
of evolution for which both VHE 7-ray and X-ray spectra of 
sufficient quality are available. Since the model is radially 
symmetric it was necessary to define a circular source area 
as an approximation to the asymmetrical Gaussian mor- 
phology fits of the published VHE 7-ray data. The radius 
was chosen such that the circle covers an area equivalent to 
the 3 a extent of the ellipse obtained from the VHE mor- 
phology fits. The values of the equivalent radii are included 
in Table [U 

Having calculated the equivalent circular extent of the VHE 



7-ray source the free parameters of the model were fixed by 
a x 2 fit to the VHE 7-ray data using the MINUIT minimiza- 
tion package ([James fc Roosl . I1975D . Note that only sta- 
tistical errors (1 a) of the VHE 7-ray data are taken into 
account. The optimized parameters and the predictions of 
the VHE 7-ray fluxes are presented in Tables [5] and [3] The 
modeled SEDs are confronted with observational X-ray, 7- 
ray (where available) and VHE 7-ray data in Fig. |3l In the 
VHE 7-ray range spectral points were available, while in 
X-rays and 7-rays the shown bands correspond to the pub- 
lished power-law fit with the parameters (index and nor- 
malization or energy flux) assumed to be uncorrelatcd. The 
confidence bands of the fitted model were calculated using 
Gaussian error propagation and take the errors of the fit 
parameters as well as their correlations into account. The 
level of the quoted error bands is adapted individually for 
the different sources for reasons of enhanced visibility (see 
caption of Fig. [3]). 

The comparison shows that the 7-ray and VHE 7-ray 
data are reasonably well described by the model (violet 
bands) , whereas the X-ray data are strongly overestimated. 
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However, this is expected since the extraction region for the 
spectrum determination is usually much smaller in X-rays 
than in VHE 7-rays. The resulting mismatch between the 
lepton population used for the modeling and the one ob- 
served in X-rays has to be taken into account. In the fol- 
lowing, we use the relation between the spatial extent of the 
lepton population and their age, given in Ea. llll Starting 
from the given size of the X-ray spectrum extraction region 
we determine the corresponding maximum age ii Gp t.,max of 
the leptons producing the emission. In the next step, we 
re-calculate the amount and energy distribution of the lep- 
tons contained within this regioro by using T — <i op t.,max 
as a starting point for the evaluation of Eq.[Tni The result- 
ing modified SEDs with their corresponding uncertainties, 
shown as blue dotted bands in Fig. [3J are clearly in bet- 
ter agreement with the observational data. A quantitative 
comparison between modeled and measured values of the 
X-ray flux is given in Table H] 

Our study indicates that, with the above introduced adap- 
tion of the spatial extent of the lepton population respon- 
sible for the observed X-ray emission, the presented time- 
dependent model allows to reproduce the broad-band non- 
thermal emission of PWNe of different evolutionary states. 
Hence, transferring the model to yet unidentified VHE 7- 
ray sources which are considered candidates for middle- 
aged PWNe, the predictions of the X-ray flux can help 
to evaluate whether deep X-ray observations with current 
satellites may be beneficial for a firm identification of the 
nature of such sources. 

4. Conclusions 

Motivated by the large number of yet unidentified VHE 7- 
ray sources suspected to be evolved PWNe, we developed 
a time-dependent lcptonic model suitable to calculate the 
non-thermal emission from PWNe of different ages. The 
presented model allows to study the expected photon SEDs 
evolving with the age of the PWN. Our study yields addi- 
tional support for the conception that evolved PWNe are 
still bright in VHE 7-rays, while their X-ray emission is 
largely suppressed and hence difficult to detect. 
Moreover, in this work we investigate the contribution of 
leptons of different epochs to the current photon SED in 
detail. In particular for older PWNe it is necessary to 
take into account the time dependence of lepton injec- 
tion and cooling effects in order to explain the observed 
VHE 7-ray emission, whereas the X-ray emission (espe- 
cially in the vicinity of the pulsar) is dominated by the 
young lepton population. 

Finally, we tested whether our model can be used to predict 
the X-ray flux of an unidentified source based on the VHE 
7-ray detection. We selected three representative PWNe of 
different evolutionary states and fixed the free parameters 
of the model by a fit to observational VHE 7-ray data. 
The comparison of modeling results and observational data 
shows that the emission in the high-energy range can be 
predicted successfully. 
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2 In a full three-dimensional treatment it would be necessary 
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